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Abstract. 

HETE-2 has provided new evidence that gamma-ray bursts may evolve with redshift Hi. We 
investigate the consequences of this possibility for the unified jet model of XRFs and GRBs [5\. We 
find that burst evolution with redshift can be naturally explained within the unified jet model, and 
the resulting model provides excellent agreement with existing HETE-2 and BeppoSAX data sets. 
In addition, this evolution model produces reasonable fits to the BATSE peak photon number flux 
distribution - something that cannot be easily done without redshift evolution. 



INTRODUCTION 

Most objects at cosmological distances (stars, galaxies, AGN) display evolution of their 
observable properties with redshift. Since gamma-ray bursts (GRBs) are thought to 
originate in the core-collapse of massive stars and are observed over a wide range in 
redshift, it is reasonable to suppose that they might also evolve. 

Since the advent of rapid GRB localizations with BeppoSAX and HETE-2, and the 
consequent follow-up observations, over 30 GRBs have reported redshift measurements. 
Using 9 BeppoSAX bursts from [1], plus an additional 11 bursts localized with HETE- 
2, ^ have been able to strengthen earlier indications that GRBs are brighter at higher 
redshifts. 

A uniform-jet model proposed by f3\ has been shown to provide a unified picture 
of GRBs and XRFs. Considering all bursts to have a "standard energy" but a range of 
jet opening solid-angles that spans five orders of magnitude, this unified jet model can 
account for the full range of observed burst properties seen by HETE-2 . Here we extend 
this unified jet model to account for redshift evolution and show that it can also explain 
the observed properties of BATSE bursts. 



OBSERVATIONS OF GRB EVOLUTION 

Analysis of the BATSE catalog has revealed evidence that GRBs may evolve strongly 
with redshift. The use of redshift estimators based on burst variability has provided 
evidence that bursts are intrinsically brighter at larger redshifts than at smaller 0,0]. 
Analyzing a set of 9 GRBs with spectroscopically determined redshifts observed by the 
BeppoSAX satellite, [IJ also claim evidence for an increase in the isotropic-equivalent 




FIGURE 1. Left: Cumulative distribution of Ey values reported by 1 2], shown with the best-fit Gaussian. 
Right: Displacement of Ey values from the best-fit central value as a function of redshift. From ||5|. 



energy (^iso) with redshift, but they did not include a discussion of possible threshold 
selection effects. 

Recent results from HETE-2 strengthen the evidence for this relationship. Figure 2 
from shows the isotropic-equivalent energies ^iso and luminosities Liso as a function 
of redshift for the HETE-2 and BeppoSAX events. After correcting for threshold effects, 
iM find that ^iso is correlated with redshift at the 5.1% confidence level, and Liso is 
correlated with redshift at the 0.9% confidence level. The observed relationship goes 
roughly as E^^o ~ (1 

UNIFIED JET MODEL SIMULATIONS 

The unified jet model of GRBs and XRFs f3\ provides a natural explanation for redshift 
evolution in GRBs. Namely, each burst exhibits a "standard energy" [Isl I2I], but the 
possible range of jet opening angles varies from fairly large values at low redshift, to 
very small values at high redshift, according to the relationship given above. That is, 
evolution in ^iso is explained by an evolution of the jet opening solid-angle, tljet- 

The burst simulations that we have implmented to test the unified jet model also 
provide a powerful way to explore models including burst evolution with redshift. For 
each burst, we obtain: (1) A redshift z by drawing from a model of the star-formation 
rate ||^, and (2) a jet-opening solid angle Qjet by drawing from specific distribution 
range in f2jet that is fixed at z = and shifts to smaller values at higher redshifts. We also 
introduce three Gaussian smearing functions to generate: (1) A spread in jet energy {Ey, 
see Figured]), (2) a spread in ^peak around the £'iso-^peak relation, and (3) a spread in the 
timescale T that converts fluence to flux. Using these five quantities, we calculate various 
rest-frame quantities (^iso, £^peakj etc.), and finally, we construct a Band function for each 
burst and transform it to the observer frame, which allows us to calculate fluences and 
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FIGURE 2. Left: Distribution of bursts in the {E'^l^, FjJ)-plane, showing the threshold curves we use 
to determine if a burst is detected by various instruments. Black points are bursts detected by the WXM, 
while gray points are not detected. Right: Distribution of bursts in the (1+z, £'iso)-plane for our best 
uniform jet with redshift evolution model. 



peak fluxes and determine if the burst would be detected by various instruments (see 
Figure |2^). 

To obtain the model presented here, we sought to roughly match the observed distri- 
bution in the (1+z, £'iso)-plane (compare Figure|2j5 with Figure 2 of |4]), which displays 
a range in Ei^o of 3 x 10"^^ to 1 .5 x 10^^ ergs at z = and 6.4 x 10^^ to 3.2 x 10^"^ ergs at 
z — 5. Assuming the faintest burst at z = corresponds to Qjet = 2;r, this translates into 
a range of Qjet of lit to 0.0125 steradians at z = and 0.0291 to 5.79 x 10^^ steradians 
atz = 5. 

The observed values of Ey are taken from and their distribution is well-fit by 
narrow Gaussian (see Figure [1^). Figure [TJ) plots the displacement of these 

values from the central value as a function of redshift and shows no evidence for 
evolution of Ey with redshift 10] ■ This rules out the possibility that redshift evolution 
might be explained by an evolution of the "standard energy". 

Since it relies on the random distribution of burst jet axes with respect to the viewing 
angle, the universal structured-jet model cannot easily accommodate redshift evolution. 
The only solution would be for the "jet energy", Ey, to evolve with redshift, but that 
would make it difficult to explain Figure 1 and the results of 0,12] • 



RESULTS 

Figure |4| compares the cumulative distributions of four observable burst quantities 
against several possible models. We find that the uniform jet model with burst evolu- 
tion can adequately describe the observed distributions of localized bursts. 

In addition, adding evolution to the uniform jet model replicates the observed distri- 
bution of peak photon number fluxes as observed by BATSE. Models without redshift 



evolution (Figure EK, dotted curve) tend to overpredict the number of high peak flux 
bursts. However, models with strong evolution (solid curve) provide excellent agreement 
with the BATSE distribution. Figure OJ? shows the differential distribution of redshifts 
(for bursts detected by WXM) for the models with and without redshift evolution. 

This model makes several predictions. Most observed XRFs are predicted to be at 
z < I. Bursts at low z have fljet > 10^^ or 0jet ~ a few degrees. We require a value of 
the "standard energy" to be ~ 5 x 10^^ ergs, or about 50 times less than the value 
reported by Thus, the fraction of Type Ic supemovae producing GRBs increases 
from ~0.1% atz = Oto ~ 10% at z ~ 5. Finally, 70% of bursts with z > 5 are detected 
by the WXM. 
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FIGURE 3. Left: Comparison of uniform jet models with redshift evolution (solid) and without (dotted) 
with the BATSE peak photon number flux cumulative distribution. Right: Predicted observed distribution 
of redshifts for the uniform jet model with redshift evolution (solid) and without (dotted). 



CONCLUSIONS 

HETE-2 has strengthened the evidence that GRBs evolve with z. The uniform jet model 
can describe XRFs and GRBs and can accommodate evolution whereas the universal jet 
model cannot. 

Confirmation of this model will require the localization of many more XRFs, the 
determination of ^peak and ^jso for many more XRFs and GRBs, and the identification 
of optical afterglows and the measurement of redshifts for these bursts. 

HETE-2 is ideally suited to localize XRFs and study their spectra, but this will be 
difficult for Swift , which has a nominal threshold of Emin ~ 15 keV and a narrow energy 
band of 15 keV < E < 150 keV. However, Swift is optimized for pinpointing X-ray and 
optical afterglows, and facilitating spectroscopic redshift measurements. Therefore, it is 
very important that the HETE-2 mission continue, even after Swift is fully operational. 
A partnership between HETE-2 and Swift can confirm or rule out GRB evolution with 
redshift. 
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FIGURE 4. Comparisons of different models with the observed cumulative distributions for Ziiso (upper 
left), Speak (upper right), 5e(2 -400 keV) (lower left) and E°^^ (lower right). The solid curve is the 
uniform jet model with redshift evolution and the dashed curve is the uniform jet model without redshift 
evolution. The dotted and dash-dotted curves are two variants of the universal or structured jet model. 
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